Comparing Simulation and Experiment of a 2D Granular Couette Shear Device 
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We present experiments along with molecular dynamics (MD) simulations of a two-dimensional 
(2D) granular material in a Couette cell undergoing slow shearing. The grains are disks confined 
between an inner, rotating wheel and a fixed outer ring. The simulation results are compared to 
experimental studies and quantitative agreement is found. Tracking the positions and orientations 
of individual particles allows us to obtain density distributions, velocity and particle rotation rate 
for the system. The key issue of this paper is to show the extent to which quantitative agreement 
between an experiment and MD simulations is possible. Besides many differences in model-details 
and the experiment, the qualitative features are nicely reproduced. We discuss the quantitative 
agreement /disagreement, give possible reasons, and outline further research perspectives. 
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I. INTRODUCTION 

The pioneering work of Reynolds in 1885 [Q and the 
more elaborate investigations by Bagnold j2| were among 
the first experiments to closely address the problem of 
granular shearing. Recently the subject of granular 
shearing has regained much interest in the physics com- 
munity due to the appearance of this process in common 
granular flows such as convection || , pipe and chute flow 
[ I, || , avalanches H 0], crack formation, and earthquakes 

In the traditional picture for shearing of a dense gran- 
ular material, grains are assumed to be relatively hard so 
that they maintain their volume and shape under applied 
forces. If shear is applied to a granular sample, in prin- 
ciple, the grains will respond elastically (i.e. reversibly) 
up to the point of failure. The response in the elastic 
regime is still an open question which is not addressed 
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because we focus on the regime of ex- 
tended deformation. Under shear, the grains will dilate 
against a normal load, up to the point of failure. Un- 
der continued shearing the system appears to approach a 
steady state, that is typically characterized by localized 
failure in narrow regions known as shear bands. An ex- 
tremely slow compaction / rearrangement can also occur 



under steady shearing |12 . However, this effect will be 
disregarded in the following. Here we are more concerned 
with the kinematics of the particles in the "quasi steady 
state" . 

Recent experiments on granular shearing have pri- 
marily focused on the force properties of the system 
H [l| 0, §]| pH, §7). Only a few experiments have 
explored the kinematics of shear zones, and these in- 
volved using either inclined or vertical chutes |j| [H| 
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or vibrated beds |19| where air flow between the particles 
may also have been important. In a single case of which 
we are aware, Buggisch and Loffelmann jl(J investigated 
the mixing mechanisms due to shearing in a 2D annular 
cell similar to the one described here. This experiment 
involved flexible boundaries, in contrast to the fixed vol- 
ume used in our work. 

The simulations presented in this study are perhaps 
unique in that they 1) match with considerable fidelity 
the parameters in a set of experiments, and 2) the sim- 
ulations and experiments yield detailed properties which 
can be mutually compared. Thus, the goals of this work 
include the deeping of insight into an important granu- 
lar system, and the opportunity to explore how well a 
class of models captures experimental observations. The 
relevance of this latter issue is underscored by a recent 
study in which diverse groups modeling flow in a hopper, 
obtained an equally diverse set ofpredictions, many of 
which did not match experiment |2l|, |2^, |2j| . 

In our experiment a 2D Couette cell filled with photoe- 
lastic polymer disks is used to study both the mean and 
statistical properties of the flow. Using particle track- 
ing techniques, the spin and transport velocity profiles 
as well as the associated density variations during steady 
state shearing can be measured. Because the particles 
are photoelastic, it is also possible to infer information 
on the local stress state of the system, a topic which is 
considered elsewhere jllj . 

The following are the key observations from these stud- 
ies: A short time after the beginning of shearing, a shear 
zone forms near the inner wall. The location of the shear 
band close to the inner wall can be attributed to the 
fact that the shear stress is highest there (decaying pro- 
portional to (r — R)~ 2 ), as was observed from previous 
simulations pi] , p5| and as is consistent with the static 
equilibrium conditions of continuum theory in cylindri- 
cal coordinates. The characteristic width of the induced 
shear band is found to be a few particle diameters, al- 
most independent of the average packing fraction. The 
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FIG. 1: (i) schematic top view of the experimental setup, (ii) 
schematic drawing of the disks close to the shearing wheel, (a) 
experimental realization of the walls, (b) realization of the 
walls in the simulation. 



mean azimuthal velocity decreases roughly exponentially 
with distance from the inner shearing wheel, and within 
the statistical fluctuations, there is shear rate invariance. 
The mean particle spin oscillates around zero as the dis- 
tance, r — R from the wheel increases, but falls rapidly to 
zero away from the shearing surface. The distributions 
for the tangential velocity and the particle spin show a 
complex shape particularly for the grain layer nearest to 
the shearing surface, indicating a complicated dynamics, 
where velocity distributions near the inner wheel are very 
wide and non-Gaussian. 



Simulations 



There are a variety of numerical studies involving 
shearing. Some of these focus on on stress-strain re- 
lations || |§ ||, |7|, || || g @, |§, and 
others deal with shear banding in specific geometries 

n, n m> n> e§ m ® m> ii- in the p resent stud y' 

we present MD simulations and investigate the kinematic 
properties of a model system which was structured so 
that its realization is as close as possible to the physical 
system discussed here, a goal partially achieved already 
in §9 . 



B. Overview of the paper 



II. EXPERIMENTAL SETUP AND 
PROCEDURE 

In this section we give enough details so that the reader 
can appreciate key aspects of the experiments and how 
the experiments relate to the simulations. The experi- 
mental setup and results are discussed in more detail in 
Refs. @, |l|, 0. The apparatus, as sketched in Fig. [l], 
consists of (A) an inner shearing wheel (with radius R), 
and (B) an outer, stationary ring confined by (C) plan- 
etary gears. In the experiments, a bimodal distribution 
of disks (D) is used, with about 400 larger disks of di- 
ameter 0.899 cm, and about 2500 smaller disks of diam- 
eter 0.742 cm. An inhomogeneous distribution is useful, 
since it limits the formation of hexagonally ordered re- 
gions over large scales, even though there might still al- 
low some short range order ]4C|| . We use the diameter, 
d = dsmaii , of the smaller disks as a characteristic length 
scale throughout this study. 

The experimental walls are fixed, corresponding to a 
constant volume boundary condition. All particles are in- 
serted into the system and the shear is applied via the 
inner wall for several rotations, before averages in the 
nominally steady state are taken. If not explicitly men- 
tioned, averages in the simulations are performed after 
about three rotations starting at t = 180 s, and extend- 
ing over three rotations, until t — 360 s. 

The mean packing fraction D (fractional area occu- 
pied by disks) is varied in the experiment over the range 
0.789 < v < 0.828. As we vary v we maintain the ratio of 
small to large grains fixed, modulo small variations due 
to the fact that particle numbers can only be adjusted 
by integer jumps. Note that the effect of the wall parti- 
cles for the calculation of the global packing fraction is 
very small. For computing the packing fraction in the 
simulations, only half the volume of the small particles 
glued to the side walls is counted, so that these boundary 
particles always contribute = 0.0047 to v. 

An important question is how the system response de- 
pends on the shearing rate, which is set by f2, the ro- 
tation rate of the inner wheel. A variation of O over 
0.0029 s -1 < < 0.09 s _1 in the experiments shows 
rate independence in the kinematic quantities, except for 
some small, apparently non-systematic variations with 
fl. A few simulations with 0.01s -1 < O < 1.0s _1 
showed clear rate independence for the slower shearing 
rates f2 < 0.1s -1 , although the situation is less clear at 
the higher end of these rates. 



The rest of this paper is constructed as follows. The 
experimental set-up is reviewed in section |ll|. The simi- 
larities and differences between the simulations and the 
physical system are discussed in section |n|. The ini- 
tial conditions and the steady state are examined in sec- 
tion [V, and the results concerning velocity- and spin- 
distributions are presented in section [v] for different den- 
sities. 



III. SIMULATION METHOD AND SIMILARITY 
TO THE EXPERIMENT 

Details of the simulations have been presented else- 
where 1 32 , and we will not repeat these. However, we 
note that the model is a soft-particle MD model. As 
noted, the parameters used in the model were chosen to 
match the experiments as reasonably as possible. Specif- 



3 



ically, the radii, static friction coefficient and density of 
the particles, and the size of the container match the ex- 
perimental values. The boundary conditions are chosen 
to mimic those in the experiment, see Sec. [H]. However, 
the "teeth" used on the inner and outer ring of the ex- 
periment are replaced by small disks with diameter d wa .n, 
see Fig. |l|. The properties of the particles and the param- 
eters for the (linear) force laws 
in Table § 



41 1 are summarized 



Property 


Values 


Diameter d sma ii, Mass m sm aii 


7.42 mm, 0.275 g 


Diameter di arge , Mass mi argc 


8.99 mm, 0.490 g 


Wall-particle diameter d wa n, 


2.50 mm 


System/disk- height h 


6 mm 


Normal spring constant k n 


352.1 N/m 


Normal viscous coefficient 7„ 


0.19 kg/s 


Tangential viscous damping jt 


0.15 kg/s 


Coulomb friction coefficient /i 


0.44 


Bottom friction coefficient fit, 


2 x 10~ 5 


Material density po 


1060 kgm" 3 



TABLE I: Microscopic material parameters of the model. 

As in the experiment, several packing fractions of the 
shear-cell are investigated in the simulations (see Ta- 
ble |n|). F° r t 00 l° w density, in the sub-critical regime, 
the particles are pushed away from the inner wall and 
lose contact, so that shearing stops. For too high den- 
sities, dilation and thus shear are hindered and the sys- 
tem becomes blocked. The intermediate regime 0.793 < 
v < 0.809 is of major interest in this study. Note that 
the range of densities that allow for the steady state 
shear flow is extremely narrow. However, we remark that 
the transition points between the three regimes quanti- 
tatively agree between experiments and simulations. 

Still, there remain some nominally modest differences 
between the experiment and the simulation, which may 
lead to differences between results for the two realiza- 
tions. The main differences are: 

• The numerical code used here accounts for a very 
weak friction with the bottom plate only, presum- 
ably smaller than reality and thus allowing less 
damping of the particles. In addition, the reduced 
friction in the experiment is achieved by powder on 
the bottom plate and this may lead to somewhat 
inhomogeneous friction between the substrate and 
particles. 

• Related to the bottom friction is a possible small 
tilt of the real particles out of plane of observa- 
tion, connected to increased tangential and fric- 
tional forces due to increased, artificial, normal 
forces. 

• The particle-wall (and also the particle-particle) 
contacts are modeled by simple linear force laws 



Global Volume 
Fraction v 



0.789 
0.791 
0.793 
0.796 
0.798 
0.800 
0.800 
0.802 
0.804 
0.805 
0.807 
0.807 
0.809 
0.810 
0.811 
0.819 
0.828 



Number 
small 



of Particles 
large 



2462 
2469 
2476 
2483 
2490 
2498 
2511 
2520 
2511 
2524 
2518 
2545 
2525 
2538 
2555 
2560 
2588 



404 
405 
406 
407 
408 
409 
400 
399 
410 
404 
412 
394 
414 
407 
399 
418 
422 



Flow Behavior 



sub-critical 



shear flow 



blocked 



TABLE II: Details of the simulation runs provided in this 
study. Mentioned are those particle numbers for which data 
were available in both experiment and simulation. The hor- 
izontal lines in the last column mark the transition between 
the sub-critical (the blocked) range of density with the shear 
flow regime. 



and thus, possibly, do not reproduce reality to the 
extent desired. More complicated non-linear or 
hysteretic or plastic models §|, HI © I3> H HI 
fl4f are far from the scope of this study. 

• In the original experiment there existed a small 
bump on the inner wheel (a deviation from the ring- 
shape, which in the end leads to a slightly larger 
effective radius of the inner wall. A larger radius 
has the strongest effect in the case of low volume 
fractions, where the particles are easily moved away 
from the inner wheel. 



There is also a difference between the way the ini- 
tial state is prepared for the experiments and the 
simulations. The starting state in the experiments 
is a nearly uniform density at the mean packing 
fraction, v. The initial state of the simulation is an 
initially dilated state, which is then compressed. 
This preparation method is described below. 



These factors apply for all the comparisons between the 
simulation- and experimental data to follow. While there 
are differences in various details, many qualitative and 
quantitative results are in agreement for the experiment 
and simulation. 
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FIG. 2: Evolution of the packing fraction v for different 
times versus radial distance (r — R) jd from the wheel in units 
of disk diameters. The upper panel shows experimental data, 
the lower one simulation results. 



IV. INITIAL CONDITIONS AND STEADY 
STATE 



In this section, we explore the initial evolution of the 
system to a nominally stationary state, characterized by 
a dilated region near the shearing wheel, with large fluc- 
tuations in local density and velocity. Before collecting 
the data in the experiment, the inner wheel ran typically 
for 60 min at the highest shearing rate, corresponding to 
at least 20 rotations of the inner wheel. In the simula- 
tion, however, the preparation had to be limited in order 
to reduce the comparatively long computation time. The 
simulations are prepared for about three rotation periods, 
because a few runs with preparation times of up to ten 
periods of rotation did not show further relaxation ef- 
fects. However, the much longer relaxation time of tens 
to hundreds of periods as used in the experiment was not 
reached, so that long time relaxation effects cannot be 
ruled out for the simulations presented here. 



extended outer ring i?p ra parc > R = 25.24 cm. While 
the outer ring is expanded, the inner ring starts to rotate 
(counterclockwise) with constant angular velocity f2 = 
0.1 s -1 . The radius of the outer ring is then reduced 
within about two seconds to reach its desired value R . 
Afterwards, the outer ring is kept fixed and the inner 
ring continues to rotate until at t = i max the simulation 
ends. Due to the constant volume in the experiment, the 
disks are inserted one by one until the desired number 
and density are reached. This difference, which can be 
seen in Fig. || affects the initial density, but should not 
influence results for the steady state. 



B. Time-evolution of density profiles 

The procedures for establishing a steady state were 
necessary because an initial, homogeneous density be- 
comes radially non-uniform as a consequence of shear- 
induced dilatancy, for both experiment and simulation. 
Starting from a fairly uniform random packing (dashed 
lines in Fig. |[ for v — 0.804), a dilated region forms close 
to the sheared inner wheel. There are minimal changes 
in the density after about 5 rotations of the inner wheel. 
Given a CPU-time of 1-2 days per rotation, we did not 
extend simulations over more than ten rotations, so that 
the true long-time behavior may not be captured here. 
Particle rearrangements have been observed over much 
longer relaxation times experimentally [fl6| . 

When making comparisions between the model and 
the experiment, it is important to keep in mind the 
following differences in obtaining local density data: 



• Fill- up procedure (Section IV A). The dashed line 
from the simulations in Fig. | is obtained after a few 
seconds of compression and shear, so that a transient 
state between the initial and the steady state of the shear 
band is visible. The experimental data are obtained 
from the static initial state, where no onset of the shear 
band could take place. 

• Averaging Procedures. The simulation data are 
averaged over full rings around the symmetry center 
of the shear-cell, whereas in the experimental system 
only radial slices that correspond to one quarter of 
the entire apparatus were observed. Even though 
averages were computed over an extended time interval, 
a systematic error due to this procedure cannot be ruled 
out. Because of possible circumferential fluctuations 
associated with this averaging process, the area under 
the experimental curves is not necessarily constant, nor 
necessarily identical to the global density. 



A. Preparation Procedure 

As noted, the preparation procedure of the simulations 
is a dilute state. Specifically, the system starts with an 



• Experimental density determination. In the experi- 
ment the local density is measured via optical intensity 
methods, where there is some uncertainty due to light 
scattering and non-linear transmission. 
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Due to these possible systematic differences between 
the local densities obtained from experiment and from 
simulation, we take the freedom to adjust the local den- 
sity data, as described below, when making comparisions 
between simulations and experiments. 



Density difference between experiment and 
simulation 



The method used to measure the local packing fraction 
in the experiment involves a calibration with some uncer- 
tainty, in addition to the fact that the real particles are 
not perfect disks as assumed in the simulation. Specifi- 
cally, data are obtained by using the fact that UV light is 
strongly attenuated on passing through the photoelastic 
disks. This technique is calibrated against packings with 
well known area fractions, such as square and hexagonal 
lattices. There are still some small systematic uncertain- 
ties in this procedure, and if one computes the packing 
fraction using the data given in the upper part of Fig- 
ure |^, a packing fraction higher than the global one is 
found. For that reason, in Fig. |^, we shift the experi- 
mental local density data downward by a constant value 
of f s hift = 0.08. 



V. CHANGING THE PACKING FRACTION 

In this section, the dependence of the local density, the 
forces, and the kinematics of the system are examined as 
a function of P, the mean packing fraction. Using this 
global density v as a parameter has led to the discovery 
of a novel transition as the system approaches a critical 
packing fraction, v c JT^| . In the experiment we found 
v c ~ 0.792 versus v c ~ 0.793 in the simulations. 

The reason for this ^-dependence is easy to understand 
by imagining what would happen if v were very low. In 
this case, grains would easily be pushed away from the 
wheel, and after some rearrangements they would remain 
at rest without further contact with the moving wall. In- 
creasing v by adding more and more grains would lead 
to the critical mean density, P c , such that there would al- 
ways be at least some grains subject to compressive and 
shear forces from the boundaries. By adding more grains, 
the system would strengthen, more force chains would 
occur, and grains would be dragged more frequently by 
the shearing wheel. If even more particles were added, 
the system would become very stiff and eventually would 
become blocked, i.e. so dense that hardly any shearing 
can take place. In the extreme limit, due to large com- 
pressive forces and deformations, permanent plastic de- 
formations might occur and brittle materials even might 
fracture. However, due to the large deformations possi- 
ble with polymeric material used in the experiment and 
due to the relatively weak forces applied, none of these 
effects is evidenced. 




10 12 
Radius «r-R)/d) 

FIG. 3: Volume fraction v, plotted against the dimensionless 
distance from the origin (r — R) / d, for different initial global 
densities v (not shifted). The open symbols give simulation 
data with v as given in the inset. The solid symbols show 
experimental data v — ^ s hift, with fahift = 0.08. 



A. Density 

We first consider the local density profiles. In Fig. |[ we 
show v vs. (r — R)/d for several v values for both exper- 
iment and simulation. The data show good quantitative 
agreement within the fluctuations between experiment 
and simulation (after the systematic shift-correction ex- 
plained above in subsection IV C ) . There is a clear differ- 
ence in density between the dynamic, dilute shear zone 
and the static outer area. From the density data, we infer 
a width of the shear zone of about 5-6 particle diameters 
- from both experiment and simulation. 



B. Velocity and spin profiles 

In this subsection, we focus on the change in the veloc- 
ity and spin profiles with changing v. In Fig. ^, we show 
data for the velocity profiles for different v from both 
the experiment and the simulation. The profiles for the 
normalized velocity, vg/(QR), show a roughly exponen- 
tial decay, although there is some clear curvature in the 
experimental data at the outer edge of the shear zone, 
where the saturation level is reached. This saturation 
level of fluctuations in the velocity is at a higher level in 
the simulations, possibly due to the systematically larger 
shear rate in simulations used to save CPU-time, or due 
to the model for bottom friction. However, the logarith- 
mic scaling over-amplifies this very small difference. 

In the experiment, the amplitude of the exponential 
term (the velocity of the particles close to the inner wall, 
vo) decays steadily to zero as v decreases towards v c . 
The simulation data show a weaker decay of the veloc- 
ity at the inner wall with decreasing density. The fact 
that vq/(£IR) <C 1 indicates that as v — > v c , either slip 



6 




0.15 
0.1 
0.05 


-0.05 
-0.1 

-0.15 
-0.2 



5 10 
Radius ((r-R)/d) 



15 



□ 












• • 






^ssi ft mmm & •§ www 






.a 


0.796 




0.811 



10 



Radius ((r-R)/d) 



FIG. 4: Velocity and spin profiles for selected packing frac- 
tions v. The solid and open symbols denote experimental and 
simulational data, respectively. 
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FIG. 5: Velocity and spin distributions close to the inner 
wall. The solid symbols denote experimental, the open sym- 
bols simulation data. 



or intermittent shear takes place at the inner wall. Only 
values of vq/(HR) = 1 would correspond to perfect shear 
in the sense that the particles are moving with the wall 
without slip. For high densities, the agreement between 
experiments and simulations is reasonable, but for low 
densities, the magnitude of the velocities differs strongly. 
This may be due to either the differences in bottom- 
or wall-friction, or due to more irregular and differently 
shaped walls in the experiments, causing more intermit- 
tency and thus reduced mean velocities. 

The experimental and the simulated profiles for the 
scaled particle spin, Sd/(£ID), evolve in a similar man- 
ner with v. Oscillations from negative to positive and 
back to negative spins are obtained, indicating at least 
partial rolling of the layers adjacent to the inner wall [ Es[ . 
The mean spins are a little higher for all the simulations 
than in the experiment, possibly due to differences in 
the bottom friction or due to differences in the shearing 
surfaces. 

The agreement in the velocity profiles is at least 
promising, given that the bottom friction may be wrong 
by a substantial amount. Especially for higher densities 
there is good quantitative agreement. Indeed, it is for 
this case that the bottom friction and wall effects are ex- 
pected to be least important, since in this regime, the 
particle-particle interaction forces are at their strongest, 
and intermittent behavior is much relatively unlikely. 
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FIG. 6: 2D probability density for v e /(Q,R) and Sd/(QD) 
for < (r — R)/d < 1. Darker regions correspond to higher 
probability densities. The upper and lower panels correspond 
to simulation and experimental data, respectively. 



C. Velocity Distributions 

^From the previous section, we infer that changing the 
packing fraction affects not only the profiles, but also 
the distributions of the velocity. In Fig. |[ we show the 
velocity distributions in a one-particle wide radial bin 
next to the inner wheel for various ^-values from the 
experiment and the simulation. 
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The data clearly show that the peaks near the origin, 
corresponding to non-rotating particles at rest, become 
weaker with increasing density. Furthermore, the regions 
with negative spin and nonzero vg grow with increasing v. 
The fact that increasing v leads to a decreasing number 
of stationary particles is not surprising. But the forma- 
tion of the second peak in the velocity distribution at 
Vg/(£IR) ~ 0.5 is not as intuitive as the small peak at 
unity. 

A key to understanding this phenomenon is contained 
in the two-variable distribution P(vg/(ClR) ) Sd/(Q,D)), 
as shown in Fig. |^ for high (right) and low (left) density v. 
The upper data are experimental and the lower data are 
from simulations. The probability is coded in grayscale 
with dark denoting higher probabilities. This figure indi- 
cates two distinct features, corresponding to two qualita- 
tively different processes. The first feature is the concen- 
tration of probability around (0,0), which corresponds 
to a state where the disks are essentially at rest, without 
spin or translation. The other feature is a concentration 
of probability around the line vg/(flR) = 1 + Sd/(ilD), 
which corresponds to non-slip motion of grains relative 
to the wheel. No-slip here means, that the particles exe- 
cute a combination of backwards rolling and translation, 
such that the wheel surface and the disk surface remain 
in continuous contact. Thus, the peak at vg/(QR) = 0, 
which is strong for low D, corresponds to particles that 
are so weakly compressed that they can easily slip with 
respect to the shearing wheel. With higher density, and 
hence greater force at the contacts between the particles 
and the shearing wheel, slipping becomes less likely and 
the combination of translation and backwards rolling is 
the preferred state. 

VI. SUMMARY 

We have reported parallel experimental studies and 
Molecular Dynamics simulations of shearing in a two- 
dimensional Couette geometry. Here, an important goal 
was to benchmark such simulations in a setting where it 
was possible to have good overlap between the parame- 
ters relevant to the simulations and the experiments. In 
most respects, the numerical results are in good qualita- 
tive, and for some quantities good quantitative agreement 
with the experimental results p6| . 

Both simulations and experiments show rate- 
independence within the statistical errors, and the 
range of rates that were studied. We have particularly 
focused on the dependence of the shearing states on 
the global packing fraction. Good agreement between 
simulation and experiment was found for the density 
profiles associated with the formation of a shear band 
next to the inner shearing wheel with a characteristic 
width of about 5 to 6 particle diameters. 

Both simulation and experiment also showed a roughly 



exponential velocity profile. However, the simulations 
did not capture the density dependence of the experimen- 
tal velocity profiles, nor some details of the shape, espe- 
cially at the outer edge of the shear band. In this regard, 
further exploration of the role played by roughness of 
the shearing surface and the effect of the particle-bottom 
friction are necessary. The former can lead to more in- 
termittent behavior, whereas the latter might explain the 
velocity-drop at the outer edge of the experimental shear 
band. 

The alternating spin profiles in experiment and simu- 
lation agreed nicely, indicating a rolling of the innermost 
particle layers (parallel to the walls) over each other. 
Outside of the shear band, rotations are not activated, 
however. From the velocity- and spin-probability densi- 
ties, a combination of rolling and sliding with the inner 
wall is evidenced. As the density decreases towards v c , in- 
creasingly more particles remain at rest - stopped by the 
bottom friction. With increasing density, more and more 
particles are dragged with the moving wall, but at the 
same time roll over each other - in layers, with strongly 
decreasing amplitude as distance from the moving wall 
increases. 



VII. CONCLUSION AND OUTLOOK 

The present study is of particular interest because of 
the intensive attempt to match as many of the detailed 
properties of the experiment as feasible by the corre- 
sponding properties in the simulation. Specifically, most 
of the parameters used in the simulation are fixed by 
experimental measurement. Nevertheless, certain prop- 
erties of the system were sufficiently complex or difficult 
to determine exactly, so that there were some differences 
between the experimental and simulational realizations. 
In this category of complex properties are friction with 
the bottom surface, and the fact that the particles were 
not perfectly uniform. In spite of these differences, all 
the features seen in the experiment were also realized 
in the simulation. In many cases, the correspondence 
between simulation and experiment were quantitatively 
correct to within a few percent. In other cases, the sim- 
ulation could be shifted appropriately so that agreement 
with the experiment was possible. Given the uncertainty 
in experimental parameters and/or small irregularities, 
this level of agreement is quite reasonable. The clear 
conclusion is that with sufficient care, MD modeling of 
a granular system can produce and predict experimental 
behavior, with the understanding that absolute quantita- 
tive agreement may be limited. Inevitably, in any exper- 
iment, small variabilities between particles or boundaries 
come into play at a sufficient level of detail. At this level 
agreement between simulation and experiment, and even 
between similar but distinct experiments must limit ab- 
solute agreement. 
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